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Abstract. In [T^, S. Montgomery-Smith provides a one dimensional model for the three 
dimensional, incompressible Navier-Stokes equations, for which he proves the blow up of 
solutions associated to a class of large initial data, while the same global existence results as 
for the Navier-Stokes equations hold for small data. In this note the model is adapted to the 
case of two and three space dimensions, with the additional feature that the divergence free 
condition is preserved. It is checked that the family of initial data constructed in [^, which 
is arbitrarily large but yet generates a global solution to the Navier-Stokes equations in three 
space dimensions, actually causes blow up for the toy model — meaning that the precise 
structure of the nonlinear term is crucial to understand the dynamics of large solutions to 
the Navier-Stokes equations. 



1. Introduction 

Consider the Navier-Stokes equations in M"^, for (i = 2 or 3, 

{dtu — Au + u ■ Vn = —Vp 
div n = 
uit=o = uo, 

where u = [u^, . . . , u'^) is the vefocity of an incompressible, viscous, homogeneous fluid evolv- 
ing in M"', and p is its pressure. Note that the divergence free condition allows to recover p 
from u through the formula 

—Ap = div (u ■ Vii). 

A formally equivalent formulation for (NS) can be obtained by applying the projector onto 

divergence free vector fields P Id — VA^^div to (NS): 

dtu -Au + F{u ■Vu) = 
u\t=o = no = P-uo- 



This system has three important features: 

(E) (the energy inequality): the L^(M^) norm of u is formally bounded for all times by that 
of the initial data; 

(I) (the incompressibility condition): the solution satisfies for all times the constraint divn = 0; 

(S) (the scaling conservation): if u is a solution associated with the data uq, then for any pos- 
itive A, the rescaled u\{t,x) An(A^t, Ax) is a solution associated with uq^x{x) *=^Auo(Ax). 
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Of course the two first properties are related, as (I) is the ingredient enabhng one to obtain (E), 
due to the special structure of the nonlinear term. 

Taking (E) into account, one can prove the existence of global, possibly non unique, finite 
energy solutions (see the fundamental work of J. Leray [Hj). On the other hand the use of (S) 
and a fixed point argument enables one to prove the existence of a unique, global solution 
if the initial data is small in scale-invariant spaces (we will call "scale-invariant space" any 
Banach space X satisfying ||A/(A-)||x = for all A > 0): for instance the homogenenous 

Sobolev space H2~^^ Besov spaces i?p,oo ^ for p < oo or the space BMO~^. We recall that 

def 

'"Ilp(R'') 



B" 



Li{R+;f) 



and 

II s\, 

\BMO- 



def 



supft^||e*^/||ioo + sup R-U [ |e*^/(t,y)|2(i2/)'Y 

i>0 V xeR'^ ^ JP{x,R) ' J 

R>0 



where P{x, R) = [0, R"^] x B{x, R) and B{x, R) denotes the ball of of center x and radius R. 
We refer respectively to [6],[T] and [10] for proofs of the wellposedness of (NS) for small data 
in those spaces. When d = 2, the smallness condition may be removed: that has been known 
since the work of J. Leray (|12j) in the energy space (which is scale invariant in two space 
dimensions), and was proved in [7], [8] for larger spaces, provided they are completions of the 
Schwartz class for the corresponding norm (Besov or BMO^^ norms). 

It is well known and rather easy to see that the largest scale invariant Banach space embedded 
in the space of tempered distributions is B^^. In three or more space dimensions, it is not 
known that global solutions exist for smooth data, arbitrarily large in B^^. We will not 
review here all the progress made in that direction in the past years, but merely recall a 
few of the main recent achievements concerning the possibility of blow up of large solutions. 
Recently, D. Li and Ya. Sinai were able in |13J to prove the blow up in finite time of solutions 
to the Navier-Stokes equations for complex initial data. We note that, as for the system that 
we construct in the present paper, the complex Navier-Stokes system does not satisfy any 
energy inequality. Before that, some numerical evidence was suggested to support the idea 
of finite time blow up of (NS) (see for instance [T5] or [9]). On the other hand in [2J a class 
of large initial data was constructed, giving rise to a global, unique solution; this family will 
be presented below. Another type of example was provided in [3]. It should be noted that 
in both those examples, the special structure of the equation is crucial to obtain the global 
wellposedness. In |14j . S. Montgomery-Smith suggested a model for (NS), with the same 
scale invariance and for which the same global wellposedness results hold for small data. The 
interesting feature of the model is that it is possible (see y[4j) to prove the blow up in finite 
time of some solutions. The model is the following: 

(TNSi) dtu - Au = {u'^) inM+xM. 

The main ingredient of the proof of the existence of blowing-up solutions consists in noticing 
that if the initial data has a positive Fourier transform, then that positivity is preserved for 
the solution at all further times. One can then use the Duhamel formulation of the solution 
and deduce a lower bound for the Fourier transform that blows up in finite time. We will not 
write more details here as we will be reproducing that computation in Section [2j 
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In this paper we adapt the construction of [14J to higher space dimensions. In order to have a 
proper model in higher dimensions it is important to preserve as many features of the Navier- 
Stokes equations as possible. Here we will seek to preserve scaling (S) as well as the divergence 
free condition (I) (as we will see, condition (E) cannot be preserved in our model). This 
amounts to transforming the nonlinear term proposed in [14] (see Equation (TNSi) above) in 
such a way as to preserve both the positivity conservation property in Fourier space and the 
incompressibility condition. This is in fact a technicality which may be handled by explicit 
computations in Fourier space; actually the more interesting aspect of the result we obtain is 
that the initial data constructed in [2] to show the possibility of global solutions associated 
with arbitrarily large initial data actually generates a blow-up solution for (TNS3). This, joint 
to the fact that we are also able to obtain blowing-up solutions in the two dimensional case, 
indicates that proving a global existence result for arbitrarily large data for (NS) requires 
using the energy estimate, or the specific structure of the nonlinear term - two properties 
which are discarded in our model. 

Let us state the result proved in this paper. 

Theorem 1. Let the dimension d he equal to 2 or 3. There is a bilinear operator Q, which 
is a d-dimensional matrix of Fourier multipliers of order one, such as the equation 

{dtu - Au = Q{u, u) in M+ x 
divu = 
n|t=o = uo 

satisfies properties (I) and (S), and such that there is a global, unique solution if the data is 
small enough in BMO~^. Moreover there is a family of smooth initial data uq, which may be 
chosen arbitrarily large in B^^, such that the associate solution of (TNSd) blows up in all 
Besov norms, whereas the associate solution of (NS) exists globally in time. 

The proof of the theorem is given in the sections below. In Section [2] we deal with the two 
dimensional case, while the three dimensional case is treated in Section [3l in both cases we 
present an alternative to the bilinear term of (NS), which preserves scaling and the divergence 
free property, while giving rise to solutions blowing up in finite time, for some classes of initial 
data. The fact that some of those initial data in fact generate a global solution for the three 
dimensional Navier-Stokes equations is addressed in Section [H 

Remark 1.1. We note that the method of the proof allows to construct blowing up solutions 
for the hyper-viscous case, meaning for equations of the form 

( dtu - A"n = Q(n, u) in M+ x R*^ 
< divu = 

[ U\t=o = Uo 

where a > I and A"/ = ^~^(|CP"/(0)- Indeed, the only important feature in order to 
construct blowing-up solutions by the method of [14] is that the system written in the Fourier 
variable, preserves the positivity of the symbol, and so, the positivity of u^{t,£,) if Uq{(,) > 0, 
for any j G {l,...,d}. 

Remark 1.2. In the two dimensional case, it might seem more natural to work on the vorticity 
formulation of the equation: in 2D it is well known that the vorticity satisfies a transport- 
diffusion equation, which provides easily the existence of global solutions for any sufficiently 
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smooth initial data. An example where the vorticity equation is modified (rather than (NS)) 
is provided at the end of Section [H below. 

2. Proof of the theorem in the two-dimensional case 

In this section we shall construct the quadratic form Q, as given in the statement of Theorem[Tl 
which allows to construct blowing up solutions for the (TNS2) system. 

Let us consider a system of the following form: 

J dtu — An = Q(n, u) — Vp 
\^ divu = 0. 

Taking the Leray projection of this equation, we obtain 

dtu- Au = FQ{u,u). 

We wish to follow the idea of the proof of [T^ , thus to construct Q = P Q as a matrix 
of Fourier multipliers of order 1, such that the product PQ preserves the positivity of the 
Fourier transform. We define Q{u,u) as the vector whose j-component is, for j E {1,2}, 

(2.1) (Q(n,n))^=^g„(Z))(nV), 

i 

and we impose that qij{D) are Fourier multipliers of order 1. For example, let us simply 
choose 

Q(e) = ieii66<o( 1 I )• 

Recalling that 




we easily obtain 

so all the elements of this matrix are positive. 
The Duhamel formulation of (TNS2) reads 

where we have denoted by Qiji^,) the matrix elements of P Q{(,). 

It is not difficult to see that all the usual results on the Cauchy problem for the Navier- 
Stokes equations hold for this system (namely results of f6],[T] and \T0\ as recalled in the 
introduction). Moreover it is clear that if the Fourier transform of uq is positive, then that 
positivity property holds for all times. 

Now let us construct a data generating a solution blowing up in finite time. We will be 
following closely the argument of [H], and we refer to that article for all the computational 
details. We start by choosing the initial data uq = (mq, Uq) such that Uq > 0, and the support 
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of Uq lies in the second and fourth sector of the complex plane, that is the zone where ^1^2 < 0; 
we also suppose this spectrum is symmetric with respect to zero (and to fix notation, that the 
support of Ug intersects the set > 1/2, for j G {1, 2}). Taking into account the divergence 

free condition which states that u^i^,) = —^^^u^iO^ deduce that U2 is supported in the 

same region as ui and is also nonnegative. 

Let us denote by A the norm of uq (which will be assumed to be large enough at the end), 
and let us write uq = Awq. The idea, as in [14J, is to prove that for any A; G N and j G {1, 2}, 

(2.2) u^{t,0 > A''e-'''2'-^(''~'h,^,,w',^\0 

where we have written, w^'-^ = {w^'-^y'' and is the restriction of wq l|g^,|>i/2 to the second 
sector of the plane. Finally the time tj. is chosen so that Iq = and tf^ — > 2^"^^ log 2. 
Notice that lim tk = log2^/^. The result (j2.2p is proved by induction. Suppose that (j2.2p is 

fe^oo 

true for k — 1 (it is clearly true for k = 0). Due to the positivity of uq, we can write 

ft 

> / e-(*"^)l«l'gjj(e)(2^'(s,0*2^'(s,0)f^s 
Jo 

and using the induction assumption, along with the support restriction of il'q"^'-' , we find that 
n^-(t,0 > re-(*-^)l«l'g,,,(0(^''-'«fc_i(.))2d.u5o'"'''*^^o^'''(0 

JO 

Jo 

where ak{t) = 2^~'^^'^''~^^lt>t^.. But Wq^^'-' * Wq~^'-' = Wq'-' , and on the support of Wq'-' we 
have QjjiO — C2^. The induction then follows exactly as in |14j . 

Once (|2.2p is obtained, the blow up of all ^ norms follows directly, noticing that W'^(too) 
can be bounded from below in ^ by C(yle~*°°2~^)2'°2(*+^)'^, which goes to infinity with k 
as soon as Ae~^°°2~^ > 1. That lower bound is simply due to the fact that (calling the 
usual Littlewood-Paley truncation operator entering in the definition of Besov norms) 

Mtoo)\\Bs = SUp2'="||Afcu(too)||L- > SUp2'=^|Afcn(too,0)| = SUp2'="|| A^(too)||Ll 

k k k 

since A^.u(too) is nonnegative. 

Remark 2.1. One can notice that as soon as the matrix Q has been defined, the computation 
turns out to he identical to the case studied in [14j . In particular the important fact is that uq 
is nonnegative (and that its support intersects, say, the set > 1/2J. 

Remark 2.2. As explained in the introduction, it seems natural to try to improve the previ- 
ous example by perturbing the vorticity equation, since that equation is special in two space 
dimensions. Let us therefore consider the vorticity oj = diu^ — 82^} ■ As is well known, the two 
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dimensional Navier-Stokes equations can simply be written as a transport- diffusion equation 
on u>: 

dtu) + u ■ Vuj — Aid = 0, 
which can also be written, since u is divergence free, 

dtoj + di{u^u) + d2{u^uj) - Acj = 0. 

Changing the place of the derivatives, and noticing that a derivative of u has the same scaling 
as oj, a model equation for the vorticity equation is simply 

dtuj + - Aw = 0. 

This simplified model is a semilinear heat equation for which the blow-up of the solution is 
well known (see [4], [5] It is also easy to see that the argument of [14J is true for this system, 
which therefore blows up in finite time for large enough initial data with negative Fourier 
transform. One can note that the equation on u becomes 

5tn + V^A"^((curln)2) - A-u = -Vp , div'u = 0, 

which blows up but does not preserve the sign of the Fourier transform. 



3. Proof of the theorem in the three-dimensional case 

The three-dimensional situation follows the lines of the two-dimensional case studied above, 
though it is slightly more technical. The main step, as in the previous section, consists in 
finding a three-dimensional matrix Q such that the Fourier transform of the product P Q has 
positive coefficients (we recall that P denotes the projection onto divergence free vector 
fields). Let us define, similarly to the previous section, the matrix 

m = \i\Me£ 1^ 1 1 1 j , 

where £^= {^e M^, < 0, 66 < 0, |6| < min(|6M6l)} • We compute easily that 

_ / 6' + el - ei6 - 66 6' + - 66 - 66 H + H - ei6 - 66 \ 

IP Q{i) = Mee\i\-^ il + el - 66 - 66 6' + H - 66 - 66 6' + il - 66 - 66 • 

V 6' + ei - ei6 - 66 6' + il - ei6 - 66 6' + ei - ei6 - 66 / 

Let us consider the sign of the matrix elements of PQ(6- The first line of the above matrix 
is clearly made of positive scalars, due to the sign condition imposed on the components of 6 
The components of the second line may be written 

6' + 6' - 66 - 66 = il- 66 + 6(6 - 6), 

which is also positive since either ^2 and 6 are both positive, in which case 6 > 6) or they 
are both negative in which case 6 < 6- Similarly one has 

6' + ei - 66 - 66 = 6' + - 6(6 + 6), 

and either 6 > 0, 6 < 0, 6 < and 6 + 6 > 0, or 6 < 0, 6 > 0, 6 > and 6 + 6 < 0. 
So the third line is also made of positive real numbers. 
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Now that it has been checked that all coefficients are positive, we just have to follow again 
the proof of the two dimensional case to obtain the expected result, showing the blow up of 
solutions to 



where Q{u,u) = ¥Q(u,u) is the vector defined as in ()2.ip . We will not write all the details, 
which are identical to the two-dimensional case, but simply give the form of the initial data, 
which is summarized in the next proposition. 

Proposition 3.1. Let uq be a smooth, divergence free vector field such that the components 
ofuQ are even, nonnegative functions, such that the support ofuQ intersects the set \^j\ > 1/2, 
for j £ {1,2,3}. Then the unique solution to (TNS3) associated with uq blows up infinite 
time, in all Besov spaces. 

We will not detail the proof of that proposition, as it is identical to the two dimensional case 
(thus in fact to [Hj). 

Of course one must check that such initial data exists. The simplest way to construct such an 
initial data is simply to suppose it only has two nonvanishing components, say Uq and Uq, and 
that the Fourier transform of Ug is supported in 1^^^2<o while intersecting the set > 1/2. 
The divergence free condition ensures that the same properties hold for Uq (and Uq is assumed 
to vanish identically). An explicit example is provided in the next section. 

That ends the proof of the "blowing up" part of the theorem. 

Remark 3.1. Notice that in that example, the energy inequality (E) cannot be satisfied, as 
it would require that {Q{u,u)\u)i2 > 0, which cannot hold in our situation if the Fourier 
transform ofu is nonnegative. 

4. Examples of arbitrarily large initial data providing a blowing up solution 

to (tns3) and a global solution to (ns) 

In this short section, we check that the initial data provided in [2] and which allows to obtain 
large, global solutions for Navier-Stokes equations, gives rise to a solution blowing up in finite 
time for the modified three dimensional Navier-Stokes equation constructed in the previous 
section. 

More precisely we have the following result. 

Proposition 4.1. Let <p be a function in 5(M^), such that > 0, and such that (j) is even and 
has its support in the region 1^^^2<0; while intersecting the set > 1/2, for j G {1,2,3}. 
Let £ and a be given in ]0, 1[, and consider the family of initial data 



Then for e > small enough, the unique solution of (NS) associated with no,e is smooth and 
global in time, whereas the unique solution of (TNS3) associated with no,e blows up in finite 
time, in all Besov norms. 





UO,e{x) = {d2(Ps{x),-diipe{x),0) 



where 
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Remark 4.1. It is proved in [2j that such initial data has a large norm, in the sense 

that there is a constant C such that 

C'\-\oge)T^ < ||no,eb-i^ < C(-log£)i 

To prove Proposition l4.1l we notice that the initial data given in the proposition is a particular 
case of the family of initial data presented in |2j, Theorem 2, which generates a unique, global 
solution as soon as e is small enough (in [2] there is no restriction on the support of the 
Fourier transform and di(j) is simply (p). So we just have to check that the initial data fits 
with the requirements of Section [3] above, and more precisely that it satisfies the assumptions 
of Proposition 13.11 Notice that 

^-(^) = (^6^(6, £"6,6 + ^) + ^6^(6,e"6,6 - ^)) 

We need to check that Uq^ > 0, for i £ {1,2,3}, and that the Fourier support intersects the 
set > 1/2. We have 

V(6 = ^"^"yi' (-660(6, e"6, 6 ± ^),e??(6,e"6,6 ± ^),o) , 

and we have clearly the desired properties. 

This ends the proof of the proposition, and of the theorem. 
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